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Abstract 

Matrix-valued stochastic processes have been of significant importance in areas such 
as physics, engineering and mathematical finance. One of the first models studied has 
been the so-called Wishart process, which is described as the solution of a stochastic 
differential equation in the space of matrices. In this paper we analyze natural exten¬ 
sions of this model, and prove the existence and uniqueness of the solution. We do 
this by carrying out a Picard iteration technique in the space of symmetric matrices. 
This approach takes into account the operator character of the matrices, which helps 
to corroborate how the Lipchitz conditions also arise naturally in this context. 
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1 Introduction 

Bru [2] introduced the so called Wishart process, which is specified by the following stochastic 
differential equation (SDE) valued in the space of symmetric dxd matrices for certain values 
a in the so-called Wallach set (i.e. a G {1, 2,..., d—1} [J[d—1, oo)) and some initial condition 
2 ^ 0 , 

dXt = + dBj\fXt + aldt, 

where B^ is & dxd matrix with each entry being a Brownian motion, all of them independent. 
If Xt belongs to the set of positive semidefinite symmetric matrices, one can properly define 
y/jQ for each t > 0, as y/JQ = Uty/A^U'f, where UtAtUf is the spectral decomposition of Xt, 
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and \/A^ represents the diagonal matrix where the diagonal is given by the sqnare roots of 
the eigenvalnes increasingly ordered. 

The work in [2] has motivated several stndies and applications, let ns mention for example 
[ini[l3l[l2]. Under certain conditions, Brn tells ns that snch an eqnation has a nniqne strong 
solntion. To argne that, Brn |2] appeals hrst to the fact, taken from [15], that the sqnare root 
is an analytic fnnction in the space of symmetric positive matrices. Then, she refers to the 
resnlt in Ikeda and Watanabe [8] abont the nniqne existence of a solntion of a vector-valned 
diffnsion to conclnde the nniqneness and existence of the solntion. Here we propose taking 
a different ronte, where one needs some resnlts on matrices. 

A more general model is given by 

dXt = g{Xt)dBJ{Xt) + f{Xt)dBfg{Xt) + b{Xt)dt, ( 1 ) 

where g, f and b are matrix valned fnnctions acting on matrices. One example is when one 
takes M —)■ M fnnctions and nses spectral decomposition to obtain matrix-valned fnnctions; 
actnally, we were motivated to stndy eqnation ([T]) after seeing this point of view in [5]. 

In this paper we propose nsing the Picard iteration method to stablish the existence of the 
solntion, as well as exploiting the operator character of the matrix to stndy the eqnation. 
We think that dealing directly (instead of passing throngh vector-valned diffnsions) with 
the matrix-eqnation can be nsefnl to obtain more insight into matrix diffnsions. Thus, 
we hope that our results help to complement theory already developed in papers such as 
[DEI El HI El El ESI El [16]. To carry out the proof, we develop a few results suited to handle 
matrix-equations and which are of independent interest. 

One important issue that we are not studying here is the so-called time of collision] one 
might read more about this in [ID E] El E] • 


2 Preliminaries 

Let Sdxd be the set of symmetric matrices, and the positive semidefinite ones. Let Bt 
he a. d X d Brownian motion (i.e. a matrix filled with independent Brownian motions). We 
will focus on the following SDE valued in Sdxd- 

dXt = g{Xt)dBJ{X,) + f{X^)dB^g{X^) + b{Xt)dt, 

with initial condition Xq G and where g, f and b are Sdxd Sdxd functions. However, 
we will be more interested in considering M —)■ M functions to construct a diffusion. In this 
case the following consideration is taken for R —)■ M functions g, f and b. If A G Sdxd, by 
g{A) we mean Hg{A)H^, where HAH^ is the spectral decomposition of A and g{A) is the 
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diagonal matrix with the values 5 '(Ai),..., g{\d) and Ai < A 2 < ... < are the eigenvalues 
of A increasingly ordered. Under this framework it turns out that Xt is a symmetric matrix 
for all t. In [7], Chapter 6, there is a detailed study of functions acting on spaces of matrices, 
an idea which is extended in functional analysis to so-called functional calculus to dehne a 
function of an operator. 

Our aim is to investigate the condition on the functions g, f and b, under which previous 
equation has a unique strong solution. As expected, Lipschitz conditions will play a crucial 
role. Before we embark on this task, some useful results are in order. 

Definition 1 For matrices A and B, by 


A > 0 

we mean A is positive semidefinite, that is Ax > 0 for every vector x, and by 

A< B 


we mean B — A>t). 

Remark 2 The following results will be useful (we used fWf as a general reference), 
i) For symmetric matrices A, B, 


(A + Bf < 2A^ + 

ii) For a symmetric matrix A and unit vector x, 

{x^AxY < x^A^x. 

A proof of ^ can he obtained using the Cauchy-Schwarz ineguality. Indeed 
{x'^AxY = {x,AxY < {x,x){Ax,Ax) = x'^A'^Ax = x'^A^x. 

Next, we prove an analogous result of the Cauchy inequality, which will be useful. 
Proposition 3 Let {At, t > 0} in Sdxd with each entry being a continuous function. 

J Agdsx^ < tx^ J A^dsx, 


( 2 ) 

(3) 


Then 


for any unit vector x. 
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Proof. First fix a unit vector x. Now consider an equidistant partition {si, ..., s„} of [0, t] 
and set Ai := Ag. for i = 1,... ,n. Let A > 0 be the partition size and define 

n 

F{u) := 'y^^{'/AAi — u'/AI)‘^x. 

i=l 


We have 


y (\/AAi - uy/Klf = + ... + AA^ -2u (AAi + ... + AA„) W (A + ... + A) 


2=1 


-V— 

t 


Since F{u) = x'^ax — 2ux'^/3x + u^x'^tx > 0 for all u, then for the discriminant 

{—2x'^I3xY — Ax^axx^tx < 0 . 

But x'^tx = t, hence 

x^{AAl + ... + AAl)xt > {x^{AAi + ... + AA„)a;)^. 

The result follows after taking the inhnitesimal sum on both sides of the previous inequality. 


3 Existence and uniqueness 

We will use the following criterion to establish the solubility of the stochastic equation. 

Definition 4 Consider a function g : S^xd S^xd- We say that g is Lipschitz in matrix 
sense if there exists a constant c > 0 such that for any pair Ai,A 2 G Sdxd o^nd any unit 
vector X we have 


x'^ {g{Ai) - g{A2)f X < cx^{Ai - A2fx. (4) 

The following two results will be useful for Theorem [71 their proofs are left in the Ap¬ 
pendix, where it is properly dehned what we mean by the matrix stochastic integral. Notice 
that the next proposition resembles an isometry property, 

Proposition 5 Let At and Ct, t > 0 6e matrix-valued stochastic processes such that 

AgdBgCs 
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is well defined as an ltd stochastic integral. Then, for any pair of vectors x,y E 


ft \ 2 

A,^dB,^C., I X 


'0 


= E [x'^CjCsAsA^y] ds. 


Lemma 6 Let t > 0 be fixed. For continuous adapted processes At and Ct in S^xd, there 
exists /3 > 0 such that 



( E 




/ E ’ 


E 

x'^ ( 1 

AsdBsCsY / CsdBjAs) x 

< 

E 

x^ ( / AsdBsCs ) X 



\Jo 

Jo J 



\Jo J 



+ fi 

E 

x^ ( f CsdB^A^ X 




\Jo J 



for all t G [0, r]. 


Theorem 7 Suppose that g, f and b are M —)■ M bounded functions that satisfy the property 
of Definition Then, the stochastic differential eguation 

g{XfidBJ{Xfi+ f f{Xs)dBjg{X,) (5) 

Jo 

has a unigue strong solution in Sdxd- 


Xt = Xo+ / b{Xfids + 


The following proof follows the general structure of its vector analogue taken from [T5] . 
Proof. Since g, f and b are bounded, there is a constant c > 0 such that for any symmetric 
matrix A 

g{A) < cl, (6) 

and the same for / and b. 

Uniqueness. Let r > 0 be hxed and consider t G [0,r]. If Xt and U are two solutions 
of the SDE ([5]), we want to prove that (Va: G E [x'^ {Xt — Ytf x'\ = 0. 

Using inequality ([2]) we have 


E[x^ {Xt-Ytfx] = E 


x^i / {b{X,)-b{Yfi)ds + Ht{X,Y)] X 


< 2E 


x^ii {b{X,)-b{Y,))ds) +{Ht{X,Y)f]x 
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where 


Ht{X,Y):= / g{Xs)dBJ{X,)+ f{Xs)dBjgiX,)- g{Y,)dBJ{Y,)- f{Y,)dBjg{Y,) 


With Proposition [3] and the Lipschitz condition (|1]) we have that 


X 


(6(W) - 6(W))ds a: 


<TC E[x'^{Xs — YsYx]ds. 


Next we do the following for the other term, 


H,{X,Y) = Ht{X,Y)± / g{X,)dBJ{Y,)± / f{Y,)dBjg{X,) 


g{Xs)dB,if{X,)-f{Y,))+ / {f{X,)-fiY,))dBjg{X,) 


hw 


+ / (.7(W) - g{Y,))dBJ{X,) + / f{X,)dB^{g{X,) - g{Y,)) 


(7) 


H(2) 


Notice that last expression is the snm of two symmetric matrices, H^^'> and thns, to 
analyze 

E[x'^{Ht{X,Y)fx], 

we can apply ineqnality (E]) to split the previous expression into two parts, one with 
and the other with H^‘^\ After that, we can apply Lemma |6]to each part with so that 
in the end we have split it into four terms. This means that E \x'^[Ht[X, Y)Yx~\ is less than 
or equal to the sum of four terms, each one of the form 


E 


X 


g{Xs)dB,{f{Xs)- f{Ys)) 


2 

X 


Using Proposition O the Lipschitz (jl]) and the boundedness conditions ([H]), the following 
happens to each term 


X 


g{Xs)dB,{f{Xs) - f{Y,)) X 


= ['e [x^(/(W) - f{Ys)fg\X,)x] ds 

Jo 

< Cl [ E[x'^{Xs - Y,fx]ds, 

Jo 
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for some finite constant Ci. All this, together with ([7]), ends up giving that 

E [x^ (Xt - Yt)^ x]<Cr [ E [x'^ (X, - Ysf x] ds, 

Jo 

where Cr is a constant depending on r. An application of Gronwall’s Lemma hnishes this 
part of the proof, which is to say that E (X^ — Y^)^ x~\ is in fact zero for all unit vectors 

X. 

Existence. The Picard iteration technique commands us to dehne 

.Y<”' :=A'„+ /'()(.¥<“-■')*+ [‘ 

Jo Jo Jo 

and xj:^^ := Xq for all f > 0. We want to prove that there exists a stochastic process Xt 
valued in Sdxd such that 

i) xj:^'^ Xt uniformly on f G [0, r] and ii) that X* satishes the SDE ([5]). 

First, in order to prove i), with techniques already used in the Uniqueness part, i.e. 
inequality ([2]), Lemma [HI Proposition |H1 as well as the boundedness condition ([HD, we have 


E[x^{xi^^ - xiyx] 


( 8 ) 


= E 


x^ / biXi^^)ds+ / ^(Xi°))d5J(XW) + / /(XW)dSj<?(Xi°)) X 


< 2eE[x^b‘^{Xo)x] + 2^eE[x^'f{Xo)g\Xo)x] + 2^eE[x^'g\Xo)fiXo)x] < c.. 


„T p2 


T „2/ 


for all t G [0, r], where Cr is a hnite constant depending on r. 
Now, using again ([2]), we have the following inequality 






= E 

x^ (6(X(’^)) - b{Xl^-^^))ds + X 


< 2E 

x^ (6(X(”)) - X 

+ 2E 

x^ {Ht{X^^\X^^-^'^)f X 


Let us analyze the last two terms in the left hand side of the previous display. With 
Proposition [3] and the Lipschitz condition (|1]) we obtain 


E 


/ (b(xi">)-b(xl’‘-»))dA X 


< T 


E [x^(Xi’^) - X^^-^^fx] ds. (9) 
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For the iJ-term, using the same idea as for the Uniqueness part: 


= ± rg(A<"VB./(-Y<"-'') ± [‘ f(Xi’'-»)dBjg(Xi’'>) 

Jo Jo 

= f 3(Xi'‘>)dB.U(Xi'‘>) - /(Ai"-‘>)) + /‘(/(-’fi”’) - /(-Yi"-‘>))<iBj9(A'i’>>) 

Jo Jo 

+ r(9(A'<”') - 9(Ai"-'l))dB,/(A<”l) + [‘ /(A<”-‘))dBf (9 (a 1”)) - g(A<”-‘>)), 


Using inequality ([2]) and iterating the same arguments (i.e. Lemma [6], Proposition [5l Lips- 
chitz condition (|1]), boundedness condition ([6])) as in the Uniqueness part we arrive at 




< Cr 
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where we also used 
incorporating ([9]). 
Dehne now 


in the last iteration. Notice that we obtain the same inequality after 


Dn := sup 

t<T 


X 




-X 


W' 


X 


We want to prove that 

CXD 

5^F(Zl„>l/n2)<cx). (11) 

71=1 

Using the Chebyshev inequality 


n=l n=l 


So, it suffices to show that 

However, 


E(Dl) < 9 


(/?’■) 


n\ 


Dn = sup 

t<T 

< r 


x^ (&(Xi")) - b{Xi^-^^))ds + a: 

x^{b{X^^^)-b{Xi^-^^))x\ds + sup\x^Ht{X^^\X^^-^^)x\ . 


t<T 








Thus, 


Di<2 


x^((6(XW) - rfs ) +2 ( sup \x^ 


t<T 


A2 


C2 


From the Cauchy-Schwarz inequality and using the inequality ([3]), we can produce 

' |x^((6(XM)< [ Ids [ \x^{{b{Xi^^) - b{Xi^-^'^)))xf ds 

3 J Jo Jo 


<T a:'^(6(XW) - 6(Xf-i)))2xds. 


Hence, for a constant depending on r. 


ml) < 4” 


im 


n\ 


For the Fff-term, since x^Ht{X^^^ — X^"^ is a martingale, by Doob’s inequality 

E{Cl) < AE[x'^{Hr{X^^\X^'^-^'^)fx]. 

Upon the same argument as for (ITOjl . 


ml) < 4 "’ 


(Pt) 

n\ 


for some hnite constant c^\ therefore, for Cr '■= max(cl^\ 

(/Sr)- 


E{Di) < Cr 


n\ 


Since flTTll holds, by the Borel-Cantelli Lemma, 

F(liminf{Zi)„ < 1/n^}) = 1. 

n 

This says that 

for almost all a; G hi there exists N{lj) such that Dn Kl/vJ for all n > N{oj), 
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which implies that is a Cauchy sequence a.s., since 


T "v" (^) 

X X; ’x = 


n—l 


T ( v{i+i) _ j^(i) j 




i=0 


The conclusion is that there exists a process Xt with 

^(n) uniformly in [0, r], as n —)■ oo. 

It remains to prove that X^ satishes ([5]), which is achieved from the inequality 




(6(XW) - b{Xs))ds + HtiX^^\ X)j X 


<Crl3 / E[x^ - XsY x]ds, 


and taking n —)■ cxo. However, previous the inequality can be obtained repeating the same 
kind of arguments used along the proof. ■ 


4 Appendix 


In this secction we prove Proposition [5] and Lemma [6l In what follows, ||H|| represents the 
operator norm and ||A||f the Frobenius norm of a matrix A. It is well known that the norms 
in the space of matrices are equivalent because it is of hnite dimension. 

We now give some remarks regarding the matrix stochastic integral I := AgdEgCg that 
we use below. Notice hrst that I is dehned as the matrix where each {i,j) entry is given by 


d d n 

Hhj) / 

k=l r=l “'0 


As{i, k)Cs{r,j)dBs{k,r) 


where A{i,k), C{r,j) and B{k,r) are the corresponding entries of A, C and B. Therefore, 
the existence of I occurs if 


r\\A{i,k)C{r,j)\\l 


E 


{As{i, k)Cs{r,j)f ds 


uo 


< CX), 


for all i, fc, r,j G {!,..., d}. The above expression || • ||2 is a norm in a space of stochastic 
processes. 

Let us see how we can construct a sequence of matrix step processes and such 
that E[\\I^'^'> — Jp] —)■ 0, n —)■ cxo, where := A^f^dB Take precisely the two 
sequences such that 

\\A"'>(t,k)Cf">(r,j) - A(i,k)C(r,])h 0 
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as n ^ oo for all i, k, r,j G d}. By the construction of the stochastic integral in one 

dimension, we can bound to have that j) —/(i, j))^] ^ 0 as n ^ oo, which helps to 

see that — I\\p] —)■ 0. Nevertheless, by the equivalence of norms — J||^] —)■ 0. 


4.1 Proof of Proposition [5] 


i) For step processes. First of all, we can check the formula for matrix step processes. In 
this case we have 



/ r \ 2 


/n-l \2 

E 

y'^ 1 / AgdBgCg j X 

= E 

- Bs,)Cg^ X 


\Jo / 


\k =0 J 


Notice that when expanding the square and taking expectation, the cross terms are vanished, 
then we have 


E 


AgdBsCs 1 X 


n —1 




X 


k =0 


Thus, we have to analyze 

written in a compact form as 


E [a^l3( 


using the notation 


® • y (d ■ c. Cs^,Ag^, b. Cgf^x. 


Notice that fd is a matrix of independent normal r.v.s with mean 0 and variance Sk+i — Sk¬ 
it will be easy to deduce the formula by analyzing the 2-dimensional case: 


E [a^l3c/db] = E 


(oi 02 ) 


fdii fdi 2 \ f Cii C12 
fd21 /^22 / V *^21 C22 


fdn fdi2 

fd21 /d22 


bi 

b2 


= E 


(iildii + 02/^21 
Ctlfdu + 0‘2/d22 


Cii C12 A / bifdii + 62/^12 
C2I C22 / I bild21 + &2/^22 


= E[aiCubi + 0102162 + 0201261 0202262 ] (Sfc+i - Sk) 


= E 


(61 62 


Oil O 12 
021 022 


Oi 

02 


(sfc+l — Sk) 


= E[b'^ca]{sk+i - Sk) = E{x'^C'^CAA'^y){sk+i - Sk). 
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This helps to see how the formula arises for step processes. 

ii) For more general processes. 

Let A and C be matrix stochastic processes where the stochastic integral I is well dehned. 
Therefore, as mentioned above, there are approximating step processes and whose 
stochatic integral converges to / in the L 2 -norm. 

By point i) above, 

Jo 

for every n > 1. Then, we want to prove that 

\E - E [y^Px] | ^ 0, n ^ oo, (12) 


and that 

f F ds ^ [ E[x^CjCsA,A'^y]ds, oo. (13) 

Jo Jo 

For (IT^ we have 

\E[y^{I{nf-P)x]\ = E[\y^{{I{n)-I)I{n) + I{I{n)-I))x\] 

< E [|/(/(n) - I)I{n)x\] + E [|//(/(n) - /)a;|] 

< ll^lllbll(i^^[||/(n)-/||||/(n)||]+i?[||J||||J(n)-/^^ 

< ii^iiiii/ii {vmnnw] + VmE) v^[ii w - in 

Hence, we obtain flT^ . because i?[||/(n) — J|p] —)■ 0 as n —)■ cx). 

For ([13]), we need to calculate 

[‘ E [x^ - Cf^C.AX) »] ds. 

Jo 

Observe that we need to calculate 


E[ f «(s)a 2 (s)ci(s)c 2 (s) - ai(s)a 2 (s)ci(s)c 2 (s))ds], 

Jo 

where ai(s) and 03 ( 5 ) are arbitrary entries of Ai^\ c”(s) and c^is) of and similarly 
without the the index n, i.e. ai(s) represents an entry of Ag. 
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After adding and substracting 01 ( 5 ) 02 ( 5 ) 01 ( 5 ) 02 ( 5 ) we can split into two terms. Let us 
elaborate one of them, the other one is similar. We have that 


(01(5)02(5)01(5)02(5) - Oi( 5 )o 2 ( 5 ) 0 i( 5 ) 02 ( 5 ))d 5 


uo 
= E 


(01(5)01(5) — Oi(5)0i(5))o2(5)02(5)(i5 


Loo 


< E 


{ai{s)ci{s) — ai{s)ci{s)yds\ / (02 (5)02 ( 5 )) 2 (i 5 


r 

r y 1 


r y ] 

h 

/ ( 01 ( 5)01 ( 5 ) — Oi(5)ci(5))^(i5 

Jo 

E 

/ (a^(5)0^(5))M5 

Jo 


where we used the Cauchy-Schwarz inequality twice, one for the integral and another one 
for the expectation. Since i?[/J(oi( 5 ) 01 ( 5 ) — oi(5)oi(5))^d5] vanishes as n —)■ 00 , we obtain 
03. and therefore the result. 

4.2 Proof of Lemma [6] 

Dehne Mt := AsdBsCg- Since Mt + Mj is symmetric, {Mt + Mjy is positive semidehnite, 
that is 

0 < (Ml + Mjf = + {Mjf + MtM^ + MjMt, 

then 

-{MtMj + MjMt) < + (M^y. 

So that 

-E [x^^MtMj + MjMyx] < E[x^M^x\ + Elx^iM^fx] 

< \E[x'^M^x]\ + \E[x'^{Myyx]\ . 

Now, for each f G [0, r], we can hnd at > 0 such that 

at \E [x^{MtMy + MyMt)x] \ < \E[x^M^x] \ + \E[x'^{Myyx] \ . 
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From the continuous trajectories of Bt, we have that a : [0,r] —)■ (0, cx)) is actually continu¬ 
ous. Let us then define 5 := min^gjo^T-] cto Then 

E + Mjfx] = E[x^Mfx\ + E[x'^{M^fx] 

+E[x^{MtM^ + M^Mt)x] 

< \E[x^M^x] I + \E[x^{M^fx] I 

+1 {lElx'^M}^ + \Elx^(M^y-x]\} . 

Dehning /9 := 1 -|- 5“^ gives the inequality. 


References 

[1] Bru, M.-F. (1989). Diffusions of Perturbed Principal Component Analysis. Journal of 
Multivariate Analysis 29, pp. 127-136. 

[2] Bru, M.-F. (1991). Wishart Processes. Journal of Theoretical Probability 4(4), pp. 725- 
751. 

[3] Cuchiero, C., Filipovic, D., Mayerhofer, E., Teichmann, J. (2011). Affine processes on 
positive semidehnite matrices. The Annals of Applied Probability 21(2), pp. 397-463. 

[4] Donati-Martin, C., Doumerc, Y., Matsumoto, H. and Yor, M.(2004). Some properties of 
the Wishart processes and a matrix extension of the Hartman- Watson laws. Publ. Res. 
Inst. Math. Sci. 40(4), pp. 1385-1412 

[5] Graczyk, P. and Malecki, J. (2013). Multidimensional Yamada-Watanabe theorem and 
its applications to particle systems. Journal of Mathematical Physics 54(2), pp. 725-751. 

[6] Graczyk, P. and Mayerhofer, E. (2011). Stochastic analysis Methods in Wishart Theory. 
CIMPA Workshop. 

[7] Horn, R.A. and Johnson, C.R. (1991). Topics in Matrix Analysis., Cambridge University 
Press. 

[8] Ikeda, N. and Watanabe, S. (1981). Stochastic Differential Equations and Diffusion Pro¬ 
cesses, North-Holland Publishing Company. 

[9] Kang, W. and Kang, C. (2013). Transform formulae for linear functionals of affine pro¬ 
cesses and their bridges on positive semidehnite matrices. Stochastic Processes and their 
Applications 13, pp. 2419-2445. 


14 



[10] Katori, M. and Tanemura, H. (2004). Symmetry of matrix-valued stochastic processes 
and noncolliding diffusion particle systems. Journal of Mathematical Physics 45(8), pp. 
3058-3085. 

[11] McKean, H.P. (1969). Stochastic Integrals, Academic Press. 

[12] Konig, W. and O’Connell, N. (2001). Eigenvalues of the Laguerre process as non- 
colliding squared Bessel processes. Electronic Communications in Probability 6, pp. 107- 
114. 

[13] Mayerhofer, E., Pfaffel, O. and Stelzer, R. (2011). On strong solutions for positive 
dehnite jump diffusions. Stochastic Processes and their Applications 121, pp. 2072-2086. 

[14] Pfaffel, O. (2008). Wishart Processes, Ph.D. thesis. 

[15] Rogers, L.C.G. and Williams, D. (2000). Diffusions, Markov Processes and Martingales, 
Cambridge University Press. 

[16] Stelzer, R.J. (2007). Multivariate Continuous Time Stochastic Volatility Models Driven 
by a Levy Process, Ph.D. thesis. 

[17] Trujillo Rivera, E.A. (2011). El Proceso de Wishart y la Dindmica de sus Eigenvalores 
y Propiedades Distribucionales Via Cdlculo Estocdstico, M.Sc. thesis. 

[18] Tudor, C. (2002). Procesos Estocdsticos, Aportaciones Matematicas, Sociedad 
Matematica Mexicana. 

[19] Zhang, F. (2011). Matrix Theory: Basic Results and Technigues (2nd ed.). Springer. 


15 



